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ABSTRACT 

Context. Diffusion of atoms can be important during quiescent phases of stellar evolution. Particularly in the very thin inert envelopes 
of subdwarf B stars, diffusive movements will considerably change the envelope structure and the surface abundances on a short 
timescale. Also, the subdwarfs will inherit the effects of diffusion in their direct progenitors, namely giants near the tip of the red giant 
branch. This will influence the global evolution and the pulsational properties of subdwarf B stars. 

Aims. We investigate the impact of gravitational settling, thermal diffusion and concentration diffusion on the evolution and pulsations 
of subdwarf B stars. Although radiative levitation is not explicitly calculated, we evaluate its effect by approximating the resulting 
iron accumulation in the driving region. This allows us to study the excitation of the pulsation modes, albeit in a parametric fashion. 
Our diffusive stellar models are compared with models evolved without diffusion. 

Methods. We use a detailed stellar evolution code to solve simultaneously the equations of stellar structure and evolution, including the 
composition changes due to diffusion. The diffusion calculations are performed for a multicomponent fluid using diffusion coefficients 
derived from a screened Coulomb potential. We constructed subdwarf B models with a mass of 0.465 M Q from a 1 M Q and 3 M Q 
zero-age main sequence progenitor. The low mass star ignited helium in an energetic flash, while the intermediate mass star started 
helium fusion gently. For each progenitor type we computed series with and without atomic diffusion. 

Results. Atomic diffusion in red giants causes the helium core mass at the onset of helium ignition to be larger. We find an increase 
of 0.0015 M for the 1 M Q model and 0.0036 M Q for the 3 M Q model. The effects on the red giant surface abundances are small after 
the first dredge up. The evolutionary tracks of the diffusive subdwarf B models are shifted to lower surface gravities and effective 
temperatures due to outward diffusion of hydrogen. This affects both the frequencies of the excited modes and the overall frequency 
spectrum. Especially the structure and pulsations of the post-non-degenerate sdB star are drastically altered, proving that atomic 
diffusion cannot be ignored in these stars. Sinking of metals could to some extent increase the gravities and temperatures due to the 
associated decrease in the stellar opacity. However, this effect should be limited as it is counteracted by radiative levitation. 

Key words, diffusion - subdwarfs - stars: evolution - stars: oscillation - methods: numerical 



1. Introduction 

Atomic diffusion in stars can change the chemical abundances 
during different evolutionary stages. For example, the chemi- 
cally peculiar subdwarf B (sdB) stars show surface abundance 
anomalies most likely caused by diffusive processes (Fontaine 
|& Chayer|1997| l. These stars, which are believed to be low mass 
(~0.5 M Q ) core-He-burning stars, have atmospheres typically 
dominated by H and deficient in He ( |Heber|1991| >. In addition, 



Michaud et al. (2007 1 found significant differences between stel- 



lar models evolved with and without diffusion on the Red Giant 
Branch (RGB) and the Horizontal Branch (HB). Some of their 
results should also apply to sdB stars, which are commonly be- 
lieved to be Extreme Horizontal Branch (EHB) stars and the de- 
scendants of red giants. 

Some sdB stars show pulsations, and two types of pul- 
sators are distinguished; the short-period variable EC 14026 
stars (Kilkenny et al. 1997[), an d the long-period variable PG 
1716 stars ( |Green e t al. 2003). The rapid oscillations in EC 
14026 stars are interpreted as low-order p-modes, driven by 
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the ^--mechanism operating in the iron opacity bump (Charpinet 
|et al.|[T996] >. The same mechanism has been shown to excite 
long-period, high-order g-modes in models slightly cooler than 
the PG 1716 stars (Fon taine et al.|2003 i. This temperature dis- 
crepancy might be related to the stellar opacities ( Jeffery & Saio 
2006b ). The local iron enhancement necessary in the driving re- 
gion around log T « 5.3 is caused by the competing diffusion 
processes of radiative levitation and gravitational settling. 

Despite the obvious importance of diffusion in these stars, 
sdB evolutionary models did not include this process consis- 
tently up to now. Instead, it was often assumed that, as a re- 
sult of gravitational settling, the thin envelope is H-rich and the 
H-profile is arbitrarily steep. However, the envelope composi- 
tion and the precise form of the He-H transition layer are im- 
portant for the pulsation modes. Furthermore, the assumption of 
such a H-envelope is inappropriate for the newly proposed post- 
non-degenerate sdB stars by |Han et al] ( |2002| ). pu et al.| ( |2008) 
showed that these stars have a broad He-H transition layer ex- 
tending to deeper regions, where diffusion is not expected to 
work efficiently. This is in contrast to the canonical post-He-flash 
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sdB models, that inherit a steep H-profile from their progenitors 
on the RGB. 

While the inclusion of gravitational settling, thermal diffu- 
sion and concentration diffusion during stellar evolution is rela- 
tively straightforward as we show here, the implementation of 
radiative levitation is far more difficult. Still, radiative levita- 
tion has been treated consistently in some stellar evolutionary 
models, for example for Population I and II main sequence stars 
( [Richard et al.|2001||2002) and HB stars ( |Michaud et al.|2007) . 
These studies have led to important insights in the stellar struc- 
ture, evolution and abundance anomalies, see Michau d & Richer] 
( 2008 ) for a review. Although it is eventually also our aim to in- 
clude radiative accelerations in our sdB models, it is beyond the 
scope of this work. 

In order to study the mode excitation in sdB stars neverthe- 
less, different approximations have been made in the literature. 
One could simply enhance Fe uniformly through the entire stel- 
lar envelope ( |Charpinet et al.|1996| |Jeffery & S"aiol|2006a| l, but 
then the stellar evolution and pulsations would be unrealistically 
altered. A more sophisticated treatment assumes a diffusive equi- 
librium profile of Fe, see |Charpinet et al7| ( |1997] >. These authors 
argued that the equilibrium state is reached on timescales much 
shorter than the sdB evolutio nary timescale (~10 5 compared to 
~10 8 years). |Hu et al. (2008) approximated the equilibrium Fe 
profile with a Gaussian function centered around logT = 5.3. 
Such a parametric approach is convenient for the study of mode 
excitation, because the exact shape of the non-uniform Fe profile 
does not greatly affect the driving, as long as the resulting opac- 
ity bump is large enough. Secondly, Fe is only enhanced in the 
driving region, corresponding to a very thin mass-shell (<10~ 6 
Mo). Thus, such a parametric approach does not model the shape 
of the Fe distribution in detail, but is more realistic than simply 
increasing the metallicity in the entire envelope. 

In this study, we constructed sdB models with gravitational 
settling, temperature diffusion and concentration diffusion from 
the zero-age main sequence (ZAMS) to the EHB. Radiative ac- 
celerations were neglected, but in order to evaluate how our re- 
sults might be affected, we enforced Fe enhancement in one se- 
ries of calculations, in the same way as in |Hu et aL| p008). For 
convenience, we will refer to the collective effects of gravita- 
tional settling, thermal and concentration diffusion as 'atomic 
diffusion' in the rest of this paper. The work presented here 
should be considered as a necessary but intermediate step to- 
wards achieving our goal of improving the seismic models for 
sdB stars. A more detailed study, including the effects on sdB 
pulsations, will be presented in a future paper as we are currently 
implementing additional transport processes in our models (e.g., 
radiative levitation, thermohaline mixing, and turbulence). Here, 
we focus on the implementation of atomic diffusion and the qual- 
itative effects on the evolution tracks and stability of the pulsa- 
tion modes of sdB stars. 

In the next section, we introduce the method followed in 
this work. In Section [2~T| we present the diffusion equations. In 
Section 2.2 we describe the stellar evolution code and the input 
physics. In Section [3] we show calibrated solar models as a test 
for the updated diffusion and evolution code. In Section |4j we 
examine the effects of atomic diffusion in red giants, which we 
used to build diffusive sdB models presented in Section [5] The 
conclusions are discussed in Section [6] 



2. Method 

2.1. The diffusion equations 

The effects of diffusion are in general small, and this has 
led to various, sometimes unnecessary, approximations in its 
treatment, e.g. considering trace elements and binary mixtures. 
Also, the calculation of the diffusion coefficients is greatly sim- 
plified by describing the interaction between particles by a 
pure Coulomb potential with a long range cut-off. However, a 
screened Coulomb potential gives more accurate results, espe- 
cially at high densities (Paqu ette et al.|1986| . Here, we calculate 
the diffusion velocities due to gradients of pressure (i.e. gravi- 
tational settling), temperature, and concentration in a multicom- 
ponent fluid. We use the routine by Thoul et al. ( 1994 ) (hereafter 
TBL) to solve Burgers flow equations (Burgers 1969). The orig- 
inal TBL routine uses diffusion coefficients derived from a pure 
Coulomb potential with a cut-off at the Debye length. We have 
updated it to make use of more accurate diffusion coefficients 
derived from a screened Coulomb potential by Paq uette et al.| 
( fT986l l. 

For a review of the basic equations describing atomic diffu- 



sion in stars we refer to Thoul & Montalban (2007). We briefly 
recall the most important ones. The Burgers diffusion equations 
are (N equations), 
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In addition, we have two constraints: current neutrality, 

^ ZjtiiWi = 0, 

i 

and local mass conservation, 
^ mjiiiWi = 0. 



(2) 



(3) 



(4) 



In the above 2N+2 equations, p,, «,, Z,, and m, denote the par- 
tial pressure, mass density, number density, charge and mass for 
species i, respectively. N is the total number of species includ- 
ing electrons. The 2N + 2 unknown variables are the N diffusion 
velocities w,, the N heat fluxes r,, the gravitational acceleration 
g, and the electric field E. We take the diffusion coefficients, Kjj, 
Zij, z'jj and z"j, from Paquette et al. (1986jl. The system of Eqs. (fill 
- Q can be written as a matrix equation, and it can then be 
solved by LU decomposition, i.e. by decomposing the matrix in 
a lower and upper triangular matrix. The procedure we followed 
is described in detail by TBL. We note that TBL eliminated the 
concentration gradient of He by demanding charge neutrality. 
However, this is unnecessary and we find it more convenient to 
keep all ionic concentration gradients. 
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A rigorous study requires the calculation of the degree of 
ionization, and the treatment of each ion as a separate species. As 
a simplification, one often uses the mean ionic charge. However, 
we found that this gives inconsistent results in regions where 
the degree of ionization changes rapidly. Also, difficulties arise 
in the case of neutral atoms, for TBL defined the concentrations 
with respect to electrons. Furthermore, we have not included dif- 
fusion coefficients of neutral atoms, which are due to atomic 
polarizability rather than Coulomb scattering. To avoid these 
complications, we simply assume full ionization for the calcu- 
lation of the diffusion velocities, although we are aware that this 
might lead to a slight underestimation of the diffusion velocities 
( |Turcotte et al.||1998[ |Schlattl|[2002] |. For the equation of state 
the ionization of H, He, C, N, O and Ne are calculated using an 
approximate pressure ionization model ( Pols et al.| 1995 i. 

The diffusion velocity can be expressed in terms of gradients 
of ion abundances, pressure and temperature, 



Wi 
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where the coefficients dfj 9 Clp and aj follow from the TBL pro- 
cedure. Having calculated the diffusion velocities, we can now 
solve the equations for the composition changes, 
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(5) 



where /? nuc ,; is the rate of change of species i due to nuclear 
reactions, <x is the turbulent diffusion coefficient for convective 
mixing, and the last term gives the composition changes due to 
atomic diffusion. 



2.2. The stellar evolution code and input physics 

The stellar models are constructed with the stellar evolution 
code STARS develop ed by |Eggleton| (1971] [19721 [T973), and 



Eggleton et al. ( 1973 i. This code has been updated by Pols et al. 



( 1995b to make use of an equation of state that includes pres- 
sure ionization and Coulomb interaction. Stancliffe & Glebbeek 
(2008) implemented gravitational settling in STARS by consid- 
ering trace elements in a H background. Here, we solve the full 
set of Burgers equations as described in the previous section. 



Stancliffe (2006) distinguishes three different methods for a 



stellar evolution code to solve the structure equations together 
with the equations for composition changes: the fully simulta- 
neous, partially simultaneous and non-simultaneous approach. 
The STARS code follows the fully simultaneous approach, i.e. 
the composition equations are solved simultaneously with the 
structure (and the mesh-spacing). While the non-simultaneous 
and the partially simultaneous approach can encounter numeri- 
cal instabilities if the diffusion timescale becomes short, the si- 
multaneous approach is not hampered by this. The drawback, 
however, is that it is computationally expensive to follow a large 
number of nuclear species in this way. We therefore limit to the 
abundance changes of seven species: H, He, C, N, O, Ne, and 
a fictitious species gathering all the other ones. Thus, the heavy 
elements that are not explicitly followed are assumed to diffuse 
equally. 

The code uses nuclear reaction rates from the NACRE com- 
pilation (Angulo et al. 1999) supplemented by reaction rates 
from Caughlan & Fowler|(|1988|. The 14 N(p,y) 15 (D reaction rate 



is reduced to 0.64 times the NACRE value, as suggested by 
|Herwig et all ( |2006| ) and |Formicola & LUNA Collaboration] 



(2002). The neutrino loss rates are from |Itoh et aL] ( fl989] [l992]l 
The metal mixture is scaled to solar abundances ( |Grevesse &| 
|Noels|1993] ). Convection is treated with the mixing-length pre- 
scription of |Bohm-Vitense| (fT958 1. We use a mixing-length pa- 
rameter (the ratio of the mixing-length to the local pressure scale 
height) of a = l/H p = 2.1, which gives an excellent fit to the Sun 
with the input physics used here (see Secfion[3]). 

Both convective and semi-convective mixing are treated as 
diffusion processes (see Eq.|5]). It is assumed that mixing occurs 
in regions where 

V rad > V ad - <W(2.5 + 20/? + 16yS 2 ), 

where f3 is the ratio of radiation pressure to gas pressure and 
<5 0v is the overshooting parameter. Convective overshooting is 
included for stars with Mzams > 1-5 M Q , using an overshoot- 
ing parameter S OY =0.12, which corresponds to an overshooting 
length of ~Q.25H P ( |Schroder et al.|1997||Pols~et al.|1997||1998l . 
For He-burning cores, convective overshooting is also included. 
The value of <S 0V = 0.12 encompasses any additional mixing be- 
yond the formal Schwarzschild boundary, regardless of whether 
it is caused by the actual overshooting of material from the con- 
vective region or whether it is caused by other processes, e.g. ro- 
tationally induced mixing, atomic diffusion. So when these mix- 
ing processes are explicitly modelled in the code, as atomic dif- 
fusion is now, the parameter <5 ov should actually be redetermined. 
Such a study is beyond the scope of this paper, and will be per- 
formed elsewhere (Glebbeek et al., in preparation). 

Mass-loss on the RGB is described by a version of Reimers' 
law (|Reimers|1975] l based on a physical approach by Schroder 
|& Cuntz| ( |2005F 



(R/R e )(L/L e )( r eff 



wind 



(M/M Q ) 
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with the parameter r\ set to 8 x 10~ 14 M yr _I . On the EHB, mass- 
loss is ignored in order to evaluate solely the effects of diffusion. 
It is expected that the rates in that phase are below 10" 12 M G yr _1 
( |Unglaub & Bues|2001t[Vink & Cassisi|2002| . 

Due to diffusion, the H-abundance in the outer stellar lay- 
ers will exceed the initial value. If the envelope is radiative and 
mass-loss is limited, the outer layer will consist of almost pure 



H. The standard opacity tables in the STARS code ( |Pols et aL 



1995 



Eldridge & Tou t 2004, C hen & ToutpOOT) are not suit 



able here, because they only go up to X = 0.8. Furthermore, in 
the standard version of the evolution code, only the opacity ta- 
ble corresponding to the zero-age metallicity is loaded during an 
evolution run. 

Therefore, we constructed new opacity tables that are valid 
for H- and He-burning regions and for compositions up to X = 1 . 
We took the opportunity to implement the most recent values: 
OPAL opacities by Iglesias & Rogers ( 1996[) smoothly merged 
with the low -temperature opacities by Ferguson et al. ( 2005| ). 
These are then combined with updated conductive opacities 
( Cassisi et al.|2007) by reciprocal addition. The covered range 
is 2.7 < log T < 8.7 and -8 < logR = logp/7^ < 1. At the high 
end of log T and log R, where values are missing due to the non- 
rectangular form of the OPAL tables, we extrapolate linearly. 

For interpolation convenience during the evolution calcula- 
tions, the tables are made rectangular in (log T, log p/T^ , Z, XF), 
where we defined the composition variable 

XF — X — Xq — Xq. 

X is the hydrogen mass fraction, Xq and Xo are the mass frac- 
tions of the enhanced carbon and oxygen, above that included in 
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the metallicity. Thus, XF is simply X during H-burning, but fol- 
lows CO enhancement during He-burning. We note that the low- 
temperature opacities by [Ferguson et al. ( 2005 ) do not include 
enhanced CO mixtures. Thus, our tables with CO enhancement 
are valid down to log T = 4, which is fine for sdB stars. 

We obtained tables for Z = 0, 0.0001, 0.0003, 0.001, 0.002, 
0.004, 0.01, 0.02, 0.03, and 0.04. For each metallicity, we built 
tables for the compositions: {X, X c +X ) = (l-Z, 0), (0.95, 0), 
(0.9, 0), (0.8, 0), (0.7, 0), (0.5, 0), (0.35, 0), (0.2, 0), (0.1, 0), (0, 
0), (0, 0.1 + 0), (0, 0.3 + 0.1), (0, 0.4 + 0.2), (0, 0.4 + 0.4), (0, 
0. 1 + 0.9 - Z). We derived the C/O ratio from a 0.47 M sdB star 
during He-burning. Due to the sensitive temperature dependence 
of the 3a reactions, the C/O ratio is higher for more massive 
stars, but luckily, it is more or less constant within the narrow 
mass range of sdBs. Thus, we do not explicitly follow Xq and 
Xo separately (as in Eldridge & Tout|20 04 ), but we use the fact 
that the C/ O ratio is a fu nction of X c + X (as in |Pols et al.|1995| 
and |Chen & Tout|2007) . 

We also take into account the effect of metal diffusion on 
the opacity, and therefore interpolate in metallicity during the 
evolution. It should be noted that the opacity is still calculated 
for a fixed metal mixture ( |Grevesse & Noels||1993| l, except in 
the case that Fe is artificially enhanced in the sdB models. This 
is justified, since the heavy elements diffuse with roughly the 
same velocities. 

The new opacity tables are written in a format suitable for 
the opacity routine that was originally developed for the CLES 
evolution code by |Scuflaire et al.| (2008 ). This routine also calcu- 
lates accurate opacity derivatives that are necessary for the sta- 
bility analysis. 

3. Solar models 

3.1. Calibration of solar models 

An important test for a stellar evolution code is provided by the 
Sun. In the standard approach of calibrating solar models, the 
mixing length parameter a, the initial He and metal abun dance 
are adjusted to produce, at the solar age (4.57 x 10 9 yr, |Bahcall| 
et al.|l995|, mo dels with observed solar radius ( 6.9599 x 10 



IU 



cm, 



Allen 



_"1973) , luminosity (3.842 x 10 33 erg/s, |Bahcall et al. 
and surface metal fraction Z s /Z s (0.0245r )Grevesse~& 
1993[l. Because the luminosity of a stellar model is sen- 



sitive to the mean molecular weight, the He-abundance can be 
adjusted to yield the solar luminosity. The mixing length param- 
eter determines the efficiency of energy transport by convection; 
at a fixed luminosity, a smaller a leads to a larger radius and thus 
a lower effective temperature. 



Pols et al. ( 1995| l found an approximate solar model, cooler 



and fainter within 0.2% and 0.7%, for a = 2.0, Y\ = 0.2812, 
Zi = 0.0188 with older opacity tables than used here (Rogers & 



Iglesias|1992| l, and no atomic diffusion. Since the stellar radius 
is also influenced by the opacity, stellar models using different 



opacity tables could require different values of a ( |Chieffi et al. 



1995 ). Furthermore, other physical inputs, such as the inclusion 
of atomic diffusion, will influence the calibration of solar mod- 
els. Therefore, we perform a new solar calibration with the up- 



dated input physics described in Section 2.2 

We found for a = 2.1, Z, = 0.01928, Y, = 0.2725 a solar 
model with a radius and luminosity both within 0.01% of the 
Sun. The surface mass fractions at the solar age are Y s = 0.248 



and Z s /X s = 0.0245, which is consistent with Grevesse & Noels 



([1993]). The latest determination of Z s /X s = 0.0165 by |Asplund 
et al.| i |2005) , although likely to be more accurate, poses a seri- 



ous problem for helioseismology as pointed out by various au- 
thors, e.g. |Serenelli et al.| (|2004|i, IMontalban et al.| ( |2004| and 
|Christensen-Dalsgaard et al.| ( |2009[ l. Also, the opacity tables use 
the |Grevesse & Noels| ( |1993| metal mixture. We have, therefore, 
not considered the new solar abundances here. 

The same ZAMS model evolved without atomic diffusion 
has at the solar age a radius and luminosity 1.6% smaller than the 
diffusive model of the Sun. This can be understood in terms of 
the greater mean molecular weight in the core due to He-settling. 
Consequently, the nuclear burning rate is higher, giving a larger 
radius and luminosity for the diffusive solar model. 

For the remainder of this paper, we use the mixing length 
parameter a = 2.1 for all our models. However, one should 
keep in mind that a could depend on the specific physical condi- 
tions, and could therefore vary with stellar mass and evolution- 
ary phase. 

3.2. Diffusion velocities in the solar interior 

It is illustrative to compare the diffusion velocities obtained with 
the updated diffusion routine to previous results, using (i) the 
simplified formulas by |Michaud & Proffitt] ( | 1 993| > (hereafter MP) 
for diffusion velocities in a fully ionized H-He mixture, and (ii) 
the original TBL routine. We calculated the H-diffusion velocity 
throughout the solar interior using Eq. (17) of MP, where we 
used the standard solar model as described above. Fig. [T] shows 
the ratio of the MP H-diffusion velocities and our exact values. 
We evaluate the contributions due to pressure and temperature 
separately. The ratio of the concentration gradients terms is not 
shown explicitly, because this almost coincides with the ratio 
of the pressure terms. Our results are also compared to values 
obtained with the original TBL routine. It is clear that the largest 
error in both the MP and the TBL approach is caused by the 
thermal diffusion term. It has already been noted by TBL that 
this error is likely caused by their assumption of zij = 0.6. MP 
solved the Burgers equations for a H-He mixture analytically 
without the heat fluxes, and numerically with the heat fluxes. 



1.5 - 



1.4 



w T.MP / ' w T.new 
w P.TBl/ w P,new 
w T,TBl7 w T,new 



S 1.3 



1.2 - 



1.1 



0.1 



0.2 



0.3 0.4 



0.5 



0.6 



0.7 



Fig. 1. The ratio of the approximate diffusion velocities of H to 
our new exact values, as a function of the solar radius. The solid 
(MP) and long-dashed (TBL) lines give the ratio of the pressure 
terms, the short-dashed (MP) and dotted (TBL) lines indicate the 
temperature term ratios. 
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They then represented the effects of the heat fluxes by an ad 
hoc correction to the results obtained when neglecting those heat 
fluxes. It is, therefore, not surprising that their largest error is 
also in the thermal diffusion term. 

In the left panel of Fig. [2] our total H-diffusion velocity is 
compared with results by MP and TBL. Our calculations give a 
slightly lower H-diffusion velocity. This is due to the overesti- 
mation of the thermal diffusion velocity in the approximations. 
A solar model evolved with the same input parameters as before, 
i.e. a = 2.1, Z, = 0.01928, F,- = 0.2725, but using diffusion ve- 
locities from the original TBL routine, results in a radius larger 
by 0.09% and is more luminous by 0.15%. Although this effect 
is small, it becomes relevant within the desired accuracy of cali- 
brated solar models. 

Equations (18) and (19) of MP give diffusion velocities of 
trace elements in a H-He background. The deviation between 
our exact results and the MP and TBL approximations is even 
worse for heavy element diffusion, as can be seen for oxygen in 
the right panel of Fig. [2] 

4. Progenitors of subdwarf B stars 

We distinguish two types of sdB stars: the canonical post-flash 
sdB star and the newly proposed post-n on-degenerate sdB star 
( |Han et al.||2002| |Hu et al.|2007| \200Sf\ The latter type origi- 
nates from an intermediate mass progenitor (2M < Mzams ^ 
4M Q ) that ignited He quiescently, whereas the progenitor of the 
canonical sdB star is a low mass star (M ZAM s < 2 M ). To 
become an sdB star, a red giant must have lost a lot of mass 
when the core was close to He ignition. For a single star this 
requires an enhanced stellar wind (D'Cruz ~et al.||1996 l. Since 
the majority of sdBs are observed in binaries (e.g. Maxted et al. 
|2001| |Morales-Rueda et al.||2006"| ), a more natural mechanism 
for the mass loss is binary interaction, which can be either stable 



Roche lobe overflow or common-envelope ejection (Han et al. 
2002) >. Regardless of the mechanism, we just removed the de- 



sired amount of envelope to construct an sdB model. 

As representative sdB progenitors, we use a 1 M G and a 3 M Q 
stellar model that we evolved from the ZAMS to the tip of the 
RGB. At 'zero-age' the 1 M model is assumed to be chemically 
homogeneous, while in the 3 M Q model 12 C has reached equilib- 
rium through the CNO cycle. We used quasi-solar abundances, 
X = 0.70, Y = 0.28, Z = 0.02. 

For comparison we computed for each stellar mass two evo- 
lutionary series, one with atomic diffusion as described further 
in Section [2~T| and one without. In the first case, all elements are 
diffused from the ZAMS to the RGB tip. 

4. 1 . Low mass progenitor 

One of the main differences between the 1 M model with and 
without atomic diffusion is the difference in the He-core mass at 



the He-flash, a fact that was already noticed by Michaud et al. 
(2007). We find that, with atomic diffusion, the red giant has a 



can grow to a higher value before the temperature for He-ignition 
is reached. We find AM core . t i P /AT s h eii = -0. 17, a somewhat lower 
dependence than found by [Rood ( 1972] > and Sweigart & Gross 
( |1978| l, -0.23 and -0.24, respectively. This could be due to the 



fact that in our models the difference in He-abundance grows 
gradually, while Rood ( 1972] > and Sweigart & Gross ( 1978] > ex- 
amined the effect of the initial composition using evolution cal- 
culations without diffusion. Also, different input physics such as 
the nuclear reaction rates might play a role. 

In a previous study, [Hu et al. ( |2007| l found a He core mass of 
0.472 M Q at the tip of the RGB for a 1 M model. The current 
lower value is caused by the new conductive opacities by Cassisi 
et al. ( 2007 1. These authors already mention that their updated 
conductive opacities cause a different thermal stratification in the 
He core than previous results, leading to a lower He mass at the 
onset of the flash. 

We also notice small, but significant, differences in the sur- 
face composition, see Fig. [3] Even after the first dredge up some 
differences remain. At the tip of the RGB, the He surface abun- 
dance of the diffusive model is 3% lower than that of the non- 
diffusive model. The metallicity is lower by 1%. 

4.2. Intermediate mass progenitor 

For the 3 M model with diffusion, the surface abundances be- 
come unrealistic. The outer layers consist of pure H almost im- 
mediately after the ZAMS, because the envelope is radiative. In 
reality, this would be prevented by competing processes such as 
radiative levitation, turbulence, rotational mixing, mass loss etc. 
However, we are not worried by this as we will remove most 
of the envelope to construct an sdB star. In any case, after the 
first dredge up, the surface abundances become comparable to 
the model without diffusion. This is partly because diffusion has 
not been efficient in the interior, during the much shorter evolu- 
tionary timescale of the 3 M Q model . The 1 M model evolved 
from the ZAMS to the RGB tip in 1.2 x 10 10 yr, while it took 
the 3 Ms model only 3.8 x 10 s yr. Furthermore, the convective 
envelope of the 3 M Q model dredges up more nuclear processed 
material than in the case of the 1 M Q model , washing away the 
effects of diffusion at the surface. 

We find, with diffusion, M coletip = 0.4403 M G , and with- 
= 0.4367 Mo. The difference is 0.0036 M , which 



out M. 



core,tip 

is larger than for the 1 M model. At first instance this might 
be surprising, but one should realize that the reason that the 
diffusive model has a larger M core t i p is different than for the 
1 M Q model. In this case, the convective core on the MS be- 
comes larger in the presence of diffusion. This is mainly caused 
by chemical composition changes leading to an increase in the 
opacity as shown by Richard ~et al.| ( |2001| l and [Michaud et al.| 
(2004). An extended convective core can burn more H, which 
leads to a larger He core at the end of the MS, and ultimately, to 
a larger He core at the RGB tip. 



^core.tip = 0.4664 M G , while without atomic diffusion we get 5. Subdwarf B models 



^core.tip - 0.4649 M Q , where we defined the He core boundary 
to be at X = 0.10. This 0.0015 M Q difference is caused by the 
fact that, due to gravitational settling, the He-abundance at the H- 
buming shell will be lower. This implies a smaller mean molec- 
ular weight at the shell for the diffusive model. Thus, the nuclear 
burning rate will be lower, and hence the maximum temperature 
in the core increases more slowly. Consequently, the core mass 



We do not consider sdBs that are merger products. 



During the evolution up to the RGB tip, all elements were dif- 
fused if diffusion was taken into account. However, on the EHB 
we must take special care because in the absence of compet- 
ing forces gravitational settling will cause all heavy elements to 
sink. The outer layers will consist of pure H, and we know this 
is not true for sdB stars. Their atmospheres are usually He de- 
ficient, a typical value is log[iV(He)/W(H)] = -2, although it 
can vary from -4 to -1 (Sa ffer et al.|1994[|Heber & Edelmann] 
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Fig. 3. Surface abundances from the ZAMS to the tip of the RGB 
as a function of the stellar age. The dashed lines are for the model 
with atomic diffusion, and the solid lines are for the model with- 
out. The left panels are for the model with Mzams = 1 M Q and 
the right panels for Mzams = 3 M . For the diffusive Mzams = 3 
M Q model , the abundances run off the scale to the values indi- 
cated with arrows. 



2004| >. The metal abundances show a wide spread from solar 
to subsolar by a factor 100, and different metals show different 
patterns, see e.g. |Fontaine et aL (2004i and O'Toole & Heber 



(2006). Current diffusion theory, with radiative levitation and 
weak stellar winds, has difficulties explaining the observed abun- 



dance anomalies. We make no such attempt here, but we evalu- 
ate how our results are affected by certain assumptions about the 
surface abundances. We examine five cases on the EHB in order 
to disentangle the impact of different diffusion processes: 

la) No diffusion, Fe is not enhanced, 
lb) No diffusion, Fe is enhanced. 

2a) Only H and He are diffused, so the surface metallicity stays 

roughly solar (Z w 0.02). Fe is not enhanced. 
2b) Only H and He are diffused, Fe is enhanced. 
2c) All elements are diffused, Fe is not enhanced. 

We note that the differences between case la) and 1 b) were 



already examined in |Hu et ah ( 2008 1. In principle, we expect that 
case 2b) is the most realistic, because diffusion calculations in 
HB stars indicate that He settles while the metals are supported 
by radiative levitation (Mich aud et al.|1983 Michaud & Richer 
2008). It is a tentative conclusion though, since it is not clear 
to what extend the results for HB stars apply to sdB stars that 
have higher surface gravities and effective temperatures. Also, 
these authors examined HB stars with Z=0.0001, and radiative 
accelerations in a solar metallicity star are much smaller due to 
saturation of the lines. Case 2c) tells us what happens if radiative 
levitation is not effective. 



5.1. Post-flash sdB stars 

The post-flash sdB models are constructed from the 1 M Q se- 
ries. If the luminosity due to He-burning exceeds 10 4 L Q , it is 
assumed that the He-flash occurs. The post-flash model is taken 
to be a ZAHB star with the same mass, core mass and chemical 
composition as the pre-flash model. In reality, the core abun- 
dances are changed slightly due to nuclear burning during the 
flash. Detailed calculations of the flash by |Piersanti et al. (2004) 
and |Serenelli & Wei ss (2005), suggest that C increases up to 5%. 
We are not concerned about this, since the p-modes are shallow 
envelope modes and do not probe the core. However, for the g- 
modes that propagate to deeper interior regions, the effects of 
nuclear burning and convective mixing during the flash (as de- 
scribed in e.g. [Dearborn et al.|2006 and |Mocak et al.|2008 1 are 
important. We explored this in |Hu et af7| ( |2009) ) . 
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Fig. 4. Evolutionary tracks for left : the post-flash sdB star, and right: the post-non-degenerate sdB star. The lines correspond to the 
cases la) solid line, 2a), short-dashed line, 2b) long-dashed line and 2c) dotted line as described in the text. 
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Fig. 5. The effect of diffusion on the the H-profiles at different EHB ages for left: the post-flash sdB star, and right: the post-non- 
degenerate sdB star. 



For the series without diffusion, the H-envelope is removed 
just after He-ignition. The compositions are kept fixed during 
mass loss, i.e. it is assumed that the mass loss happens on a 
timescale shorter than the nuclear and diffusion timescales. The 
total amount of H left is Me = 10~ 4 M . In this way, we ob- 
tained sdB models with M core = 0.465 M Q . To obtain sdB mod- 
els with the same core mass for the series with diffusion, we had 
to remove the envelope before He-ignition. Still, the core mass 



is large enough so that He is ignited after mass-loss, see Hu et al 
J20071 . 

In the left panel of Fig. |4j we show the sdB evolutionary 
tracks in the log g - T e ff diagram for the cases la) solid line, 2a) 



short-dashed line, 2b) long-dashed line and 2c) dotted line. We 
observe the following: 

- Compared to la) no diffusion, no Fe enhancement, track 2a) 
with H-He diffusion and no Fe enhancement, is shifted to 
lower surface gravities and effective temperatures. Due to 
the outward diffusion of H, the density in the envelope de- 
creases. Consequently the envelope becomes less gravita- 
tionally bound, and the sdB star gets larger and cooler. 

- Comparing track 2a) to 2b) with H-He diffusion and with Fe 
enhancement shows that Fe enhancement has only a small 
effect on the sdB evolution as we already saw in Hu et al. 
(2008). Due to the increased opacity caused by Fe enhance- 
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ment there is a slight shift to lower gravities and tempera- 
tures. 

- Interestingly, track 2c) with all elements diffusing, and no Fe 
enhancement, is shifted to higher gravities and temperatures 
compared to la). This is because the opacity decreases as the 
heavy elements sink. Consequently the star becomes more 
compact and hotter, and apparently this effect is greater than 
what happens due to H-He diffusion in 2a). 

To illustrate the effect of time-dependent diffusion, we show 
in Fig. [5] the H-profile against the fractional mass log( 1 -M r /M t ) 
at four different ages of the sdB star: 0, 10 7 , 5 x 10 7 , and 10 8 yr 
after the ZAEHB. Without diffusion the H-profile would hardly 
change. 

After each 10 7 years of evolution on the EHB, the pulsational 
properties were calculated with the non-adiabatic stellar oscilla- 
tion code MAD (Dupret 2001 ). In the top panels of Fig. [6] we 
show the frequencies of pulsation modes < / < 2 as a function 
of the effective temperature. We plotted the results for the tracks 
with Fe enhancement in order to evaluate the mode excitation, 
namely lb), no diffusion, with Fe enhancement, and 2b), H-He 
diffusion, with Fe enhancement. We notice distinct changes in 
both the frequencies and the frequency range of excited modes. 

To understand the latter, one must realize that the driving is 
caused by an opacity bump at logT = 5.3. When H diffuses 
outwards, the sdB star is cooler and has a lower temperature gra- 
dient due to the lower density in the envelope by 



dT 
~dr 



3 KpF 
'Aa'c^' 



(6) 



It should be noted that an increase in the H-fraction tends to in- 
crease the opacity, which also affects the radiative energy trans- 
port Eq. ([6]). However, the work of |Jeffery & Saio| ( [20 06a) indi- 
cates that the effect of reducing the mean molecular weight out- 
weighs the effect of increasing the opacity, and our results here 
regarding the mode excitation support their findings. Thus the 
opacity bump will be located deeper in the star in the presence 
of H-He diffusion. For a mode to be driven, the amplitudes of the 
eigenfunctions must be significant in the driving region, which 
occurs when the last node is at a certain temperature (log T ~ 6, 
see |Hu et al.|20 08 for details). So the last node of a mode must 
also be located deeper in the star in order to get excited. This 
corresponds to modes of lower radial order and frequency. 

5.2. Post-non-degenerate sdB stars 

For this type of sdB stars, we used the 3 M Q model as progenitor. 
The envelope is removed just after He ignited quiescently in the 
core until the total mass is 0.465 M . 

In the right panel of Fig. |4] we show the sdB evolutionary 
tracks in the log g - T e ff diagram for the same cases as before. 
Here, we deduce that; 

- Compared to la), no diffusion, no Fe enhancement, track 
2a), with H-He diffusion and no Fe enhancement, is dras- 
tically shifted to lower gravities and effective temperatures. 
In the presence of diffusion the initial low H-abundance 
(X s = 0.18) in the sdB envelope cannot be maintained. 
The post-non-degenerate sdB star has a lot of H 'hidden' 
in deeper layers. When this diffuses outwards, the resulting 
shift in the logg - diagram is much more pronounced 
than for the post-flash sdB star. 

- Again we find that Fe enhancement shifts track 2b), with 
H-He diffusion and with Fe enhancement, to slightly lower 
gravities and effective temperatures compared to 2a). 
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Fig. 6. The frequencies of the stable (x) and unstable (+) modes 
with / < 2 as a function of the effective temperature. The upper 
panels are for the post-flash sdB star and the lower ones are for 
the post-non-degenerate sdB star. We used the structure models 
from evolutionary tracks lb) and 2b) for the case of no diffusion 
{left) and with diffusion (right), respectively. 



- Track 2c), with all elements diffusing, and no Fe enhance- 
ment, shows an interesting feature. It is shifted to higher 
gravities and effective temperatures compared to track 2a), 
but still has lower gravities and temperatures compared to 
la). This is in contrast to what happens in this case for the 
post-flash sdB star. It is apparent that the decrease in opac- 
ity has a smaller impact than the effect of H-He diffusion 
described in 2a). 

In the right panel of Fig. [5] one can see the effect of H-He 
diffusion on the H-profile. The abundances in the outer layers are 
changed on a short timescale. At depths of log(l -M r /M*) > -3, 
diffusion does not work efficiently on the EHB timescale (~ 10 s 
yr). The change in the H-abundance near the core at log(l - 
M r /M*) > -1, is caused by H-shell burning. 

In the bottom panels of Fig. [6] we plot the frequencies of pul- 
sation modes < I < 2 against the effective temperature for the 
cases lb) and 2b). The effects of diffusion on the pulsations are 
distinct. This is not surprising considering the great impact dif- 
fusion has on the structure of the post-non-degenerate sdB star. 
A previous study ( |Hu et al.]20 08) found that the frequency range 
of excited modes is one of the main discriminators between the 
post-flash and the post-non-degenerate sdB star. If diffusion is 
effective, this conclusion must be revised. We still see differ- 
ences in the frequency spectrum between models from different 
evolutionary channels, but a more detailed study including more 
models is necessary to quantify the differences. We plan to per- 
form such an analysis in the near future. 

We conclude that atomic diffusion causes a larger change of 
H surface abundance in the post-non-degenerate sdB star than 
in the post-flash one. So the impact of H-He diffusion on the 
evolutionary tracks and driving is larger. In contrast to the post- 
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flash case, this cannot be compensated by a possible decrease in 
the opacity by the heavy elements settling. 



6. Conclusions 

We have updated the diffusion routine by TBL to make use of 
diffusion coefficients derived from a screened Coulomb poten- 
tial ( |Paquette et al.|1986] > rather than a pure Coulomb potential. 
The improved accuracy is mainly in the thermal diffusion term, 
and this results in slightly lower diffusion velocities compared to 
TBL. Although the effect is small, it becomes important within 
the accuracy desired for calibrated solar models. We found an 
excellent fit to the Sun by evolving a 1 M with initial abun- 
dances Z = 0.01928, Y = 0.2725 and mixing length parameter 
a = 2.1 with atomic diffusion to the age of 4.57 x 10 9 yr. 

We evolved a 1 and 3 M stellar model from the ZAMS to the 
RGB tip, with and without atomic diffusion. By including diffu- 
sion, we find an increase in the He core mass at the RGB tip of 
0.0015 Mq and 0.0036 M , respectively. Surface mass fractions 
are slightly changed for the 1 M Q model, but not significantly 
for the 3 M G model. The reason is in part the shorter evolution- 
ary timescale of the higher mass star. More importantly though, 
is the more efficient mixing of nuclear processed material dur- 
ing the first dredge up. This is because for the higher mass star, 
the inward advance of the convective envelope reaches regions 
where the shrinking convective core passed through during the 
main sequence. 

From the RGB models, we constructed post-flash and post- 
non-degenerate sdB models. Radiative levitation was not in- 
cluded, but in one set of our calculations we enhanced Fe arti- 
ficially around log T — 5.3 as an approximation to the expected 
diffusive equilibrium profile (Charpine t et al.|1996 i. This allows 
us to study the excitation of the pulsation modes, while the effect 
of such an artificial Fe profile on the sdB evolution is minimal. 
Although the question remains if such models are suitable for a 
precise seismic analysis, they are very useful for a comparative 
study as presented here. 

First of all, we see significant shifts in the evolutionary tracks 
when H-He diffusion is included. If H diffuses outwards, the en- 
velope is less dense and thus less gravitationally bound. This 
results in larger radii, and therefore lower gravities and effective 
temperatures. For the post-non-degenerate sdB star the differ- 
ence is much more pronounced, because the initial (ZAEHB) 
H-abundance in the envelope was very low, X = 0.18, and H ex- 
tended to deeper layers. In the presence of H-He diffusion in our 
sdB models, the frequencies of excited modes are lower and the 
frequency spectrum is more densely spaced. Especially for the 
post-non-degenerate sdB star the effect is drastic due to the large 
structural change of the envelope. 

We also examined what happens if the metals are allowed 
to sink due to gravitational settling. The corresponding decrease 
in opacity will tend to make the star more compact and hotter, 
thus competing with the effect of He settling. In reality, how- 
ever, radiative levitation prevents the metals from sinking, so 
we expect the case of only H-He diffusion to be more realis- 
tic. It is clear that consistent modelling of radiative levitation 
is a missing piece in this study, and we intend to include this 
in future work. Still, the results presented here, are an improve- 
ment to previous sdB models that altogether neglect atomic dif- 
fusion. In particular, we find it cannot be ignored in the post- 
non-degenerate sdB stars as it leads to totally different stellar 
structures. Furthermore, for H and He, the radiative forces do 
not play a significant role compared to settling. Therefore, the 



main conclusion of our work remains valid even when other pro- 
cesses are included: The settling of He leads to a H-rich sdB en- 
velope, which leads to a lower moleculair weight. This results 
to a lower density at the same pressure and temperature, which 
changes both the evolution tracks and the stability of the modes 
in the way presented here. 
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